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Abstract 

It is shown that smooth maps / : S 3 — > S 3 contain two countable families of harmonic 
representatives in the homotopy classes of degree zero and one. 



1 Introduction 

Let / : M — > N be a smooth map between Riemannian manifolds (M, g) and (N, h). The 
energy of the map / is defined by 

where x % denote local coordinates in M and f A denote local coordinates of the point 
f(x) in N. The critical points of this energy functional are called harmonic maps [1] (in 
physical literature the model defined by (1) is called the cr-model). 

In this paper I consider the case when both M and N are three dimensional spheres 
S 3 . Let (ip, $,<£>) be standard hyperspherical coordinates on M in which ds 2 = dip 2 + 
sin 2 -?/' (di? 2 + sin 2 -^ dip 2 ), and let (\1/, G, $) be similar coordinates on the target space N . 
Assume further that the map / is 5*0(3) symmetric: 

m> = /(VO, e = 7?, $ = ^. (2) 

Then the energy functional (1) reduces to 

E(f) = 4vr r (f 2 + 2 ^-) sin 2 ^ # , (3) 

jo \ sin ip J 

and the corresponding Euler-Lagrange equation is 

(sin 2 ^/')'-sin(2/) = 0. (4) 

By the principle of symmetric criticality [2] the solutions of this equation are the critical 
points of the energy functional (1), and therefore determine harmonic maps. The reduced 
model is invariant under the discrete transformations / — > / + wn and / — ► — /. To get 
partially rid of this degeneracy I shall work in the configuration space C of equivalence 
classes defined by / ~ / iff cos / = cos /. Note that in C the two functions / and / = n — f 
are considered as distinct. 

I am interested in solutions of Eq.(4) for which the energy density 

P=f + 2^L (5) 

sin ip 



is everywhere finite. Such solutions will be referred to as regular. In C there are two trivial 
regular solutions, namely the vacua /o = and /o = tt for which the energy (3) attains 
the global minimum E — 0. 

To my knowledge the only nonconstant regular solution of Eq.(4) known in literature 
is an identity map / = ip of degree one (and the corresponding antipodal map / = 7r — ip 
of degree minus one). 

In this paper I shall prove that Eq. (4) has a countable family of regular solutions of 
degree zero and one, which will be referred to as even and odd solutions respectively. For 
even solutions 

/M=/(7T-^), (6) 

while for odd solutions 

/W0=7r-/(7r-^). (7) 

2 Numerical results 

In order to analyse solutions of Eq.(4) it is convenient to introduce a new independent 
variable x = ln(tan(^>/2)) which changes from — oo (for ijj — 0) to +oo (for ip = n), and 
define h — f — n/2. In terms of these variables Eq.(4) reads 

h" - tarih(x)h' + sin(2/i) = . (8) 

Under the reflection x — ► —x even solutions are symmetric h(x) = h(—x), while odd 
solutions are antisymmetric h(x) = —h(—x), hence it is sufficient to consider Eq.(8) for 
x > 0. The formal power series expansion near x = for odd solutions is 

h(x) = bx + 0(x 3 ), (9) 

and for even solutions 

h(x) = d + 0(x 2 ). (10) 

Since x = is a regular point of Eq.(8), these power series have nonzero radii of conver- 
gence and therefore define local solutions near x — 0. 



For regular solutions the asymptotic behaviour near x = oo is 

±h(x) = -- + ce- x + 0{e- 3x ). (11) 



I have integrated Eq.(8) numerically using a standard shooting procedure for solving 
two-point boundary value problems [3]. The idea of this method is to find such initial 
data (9) or (10) which give rise to correct asymptotic behaviour (11). Below I describe 
the results for odd solutions, i.e., the initial condition (9) (the behaviour of even solutions 
is analogous). For b > 1 the function h(x) increases monotonically, crosses h = n/2 at 
some x\ and tends to +oo (as e 2x ) for x — > oo. As b — »■ b\ — 1, the point x\ moves to 
+oo and for b = 1 one gets the first odd solution hi. Actually this solution is nothing 
else but the identity map hi = n/2 — 2arctan(e x ). This is the only solution which is 
known analytically. For b slightly below one the function h(x) after reaching a maximum 
monotonically decreases to — oo, crossing h = —tt/2 at some x^. As b — ► 63 = 0.29696... 
the point X3 moves to +00. The value 63 corresponds to the second odd solution h%. As b 
decreases to zero this behaviour repeats itself, i.e., there is a countable sequence of initial 
values b n > (n odd) for which the solution h n (x), staying within the strip (— 7r/2,7r/2), 
oscillates ^ - times around h — and goes to h n (oo) = (—l) IL 2~ir/2. 

Similarly, there is a countable set of initial values d n defining even solutions. Both for 
odd and even solutions the index n is equal to the number of zeros of h n (x). 

The numerical results are displayed in Tables 1, 2, 3, 4 and Figs 1, 2. In Tables 3 and 
4 the symbol x^. denotes the location of fcth nonnegative node of h n (x). For aesthetic 
reasons, in Figs 1 and 2, 1 used the symmetry h(x) — > —h(x) to have the same asymptotics 
for all solutions. 

3 Large n behaviour 

For large n the numerical data exhibit a striking scaling behaviour and the solutions 
seem to converge to some limiting configuration. In order to understand this remarkable 
feature it is helpful to distinguish three characteristic regions. 



Table 1: Parameters of odd solutions 



n 


b n 


Cn 


b 2 JK + 2 


C n +2/C n 


E n /(l2n 2 ) 


1 


1 


2 


11.3397 


12.1254 


0.5 


3 


0.29696052 


24.25080426 


10.7816 


10.8061 


0.65168569 


5 


0.09043935 


262.0555747 


10.7519 


10.7538 


0.66527636 


7 


0.02758136 


2818.098290 


10.7493 


10.7495 


0.66653735 


9 


0.00841249 


30293.20888 


10.7491 


10.7491 


0.66665464 



Table 2: Parameters of even solutions 



n 


d n 


Cn 


d 2 Id 2 

u nl u n+2 


c n+2/ c n 


E n /(12n 2 ) 


2 


0.33480103 


7.07540708 


11.1126 


11.0168 


0.61595784 


4 


0.10043365 


77.9483657 


10.7794 


10.7676 


0.66199596 


6 


0.03059018 


839.318687 


10.7519 


10.7508 


0.66623251 


8 


0.00932911 


9023.30712 


10.7493 


10.7492 


0.66662628 


10 


0.00284544 


96993.7060 


10.7491 


10.7491 


0.66666291 



In the inner region (I), < x< 1, h n (x) is small and therefore is well approximated 
by the solution of the linear equation 



h" - tanh(x)/i' + 2h = . 



'121 



In this region we have approximately 



h n {x) — b n gi(x) for n odd and h n (x) = d n gi(x) for n even , 



(13) 



where the functions Qi(x) (i = 1 for odd and % — 2 for even solutions) solve Eq.(12) with 
initial conditions (#i(0) = 0, g[(0) = 1) and (#2(0) = 1, g' 2 (0) = 0), respectively. As n 
goes to infinity the zeros of h n (x) tend to the zeros of the corresponding function gt(x), 
as is shown in Tables 3 and 4. 



Table 3: Location of zeros for odd solutions 



n 


X i 


x 2 


X 3 


3^4 


X 5 


1 













3 





2.230560 








5 





2.194022 


4.606835 






7 





2.190795 


4.565917 


6.982062 




9 





2.190496 


4.562291 


6.940749 


9.356924 


oo 





2.190465 


4.561920 


6.936712 


9.311532 



Table 4: Location of zeros for even solutions 



n 


X i 


x 2 


Xo 


%4 


X 5 


2 


1.035870 










4 


1.012521 


3.394669 








6 


1.010474 


3.354995 


5.770839 






8 


1.010284 


3.351483 


5.729639 


8.145811 




10 


1.010267 


3.351158 


5.725986 


8.104472 


10.520646 


oo 


1.010265 


3.351124 


5.725613 


8.100431 


10.475252 
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In the intermediate region (II), 1 <x< n, h n (x) is still small (so (13) still holds) and 
is well approximated by the solution of the linear equation 

h"-h' + 2h = 0, (14) 

hence in this region the solution is almost periodic 

h n (x) ~ A n e* sin(— x + 5i) , (15) 

where the amplitude A n is proportional to b n (or d n ) and the phase 5i depends only on 
whether the solution is odd or even. 

Finally, in the asymptotic region (III), x>n, h n (x) is no longer small but can be well 
approximated by the solution of 

h" - ti + sin(2/i) = . (16) 

The behaviour of this autonomous system determines the asymptotic behaviour of solu- 
tions of Eq.(8). Note that Eqs.(14) and (16) are translation invariant, i.e., if h(x) is a 
solution so is h(x + A). Since asymptotically ±h n (x) = — 7r/2 + c n e~ x , we deduce that in 
regions II and III approximately 

± h n (x) = s(x - ln(c n )) , (17) 

where s(x) is the solution of Eq.(16) normalized by s(x) = —it/ 2 + e~ x for x — ► oo. 
Now, consider two solutions h n+ 2 and h n . In regions II and III one has from (17) 

h n+2 (x)~-h n (x-ln(^)). (18) 

Since in region II the solution h n+ 2{x) has one zero more than h n (x), and by Eq.(15) 
the distance between two adjacent zeros is equal to 2n/\^, one obtains from (18) the 
relation 

which explains the numerical fact (see Table 1 and 2) that 

lim — = e^ w 10.7491 . (20) 
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The scaling relation for the initial parameters b n and d n may be deduced in a similar 
manner. Namely, it follows from (13) that in regions I and II one has approximately 

b d 

h n +2(x) = -^ — h n (x) for n odd and h n+ 2(x) = —\ — h n (x) for n even . (21) 
o n ct n 

Combining (21) with (15), (18) and (19) one obtains 

d l _. K _ c„ +2 _, _^ 



^/ 2 h 2 

"n+2 u n+2 



e v7 . 22 



From this it follows in turn that the amplitude of oscillations of h n (x) in region II (and 
I) behaves as A n ~ e -A. Thus, as n — >■ oo, the solution h n {x) tends pointwise to zero 
for any finite x. The limiting solution hoc = is singular (because it fails to satisfy the 
boundary condition (11) and therefore its energy density p diverges at north and south 
poles), nevertheless its energy is finite Eoq — | x 12-7T 2 . 

Actually it is easy to show that E = 8tc 2 is the upper bound for energy of regular 
solutions. To see this rewrite Eq.(8) in the form 

ti) = — sm(2h). (23) 



v cosh x J cosh x 

Multiplying this by h and integrating by parts yields an identity 



OO Art* fCO Art* 

h' 2 = -£-hM*h), (24) 



-oo cosh a; J-oo cosh a; 

which implies that on-shell 

/oo Ax 
— — (/isin(2/i) + 2cos 2 (/i)). (25) 

-oo cosh a; 

It will be shown below that for regular solutions —it/2 < h < it/2. Within this range 

the function F(h) = hsin(2h) + 2cos 2 (/i) has a maximum at h — 0, hence the maximal 

energy is 

E = 8tt r -^— = 8n 2 . (26) 

J-oo cosh a; 

Finally, using the scaling relations derived above one can also show that for large n 

Em — E„ 22L , _ s 

" ^eV7. (27) 



Eoo — E n+ 2 



Numerical fit gives the folowing empirical formula for energy 



| 2i n7V 

E n = --ce Sf , (28) 



12tt 2 3 

where c ~ 0.526 for odd solutions and c ~ 0.539 for even solutions. 



4 Stability analysis 

Now I turn to the stability analysis of harmonic maps described above. The Hessian of 
the energy functional at h is 

^( fc )(«) = 4,£(^-2cos(2/ l )^)^, (29) 

which leads to an eigenvalue problem 

cosh x cosh x cosh 3 x 

Let ££ (and resp. AJ?) denote the fcth (where k = 1,2,...) eigenfunction (eigenvalue) 
around the regular solution h n . For hi = tt/2 — 2arctan(e :r ) Eq.(30) is solved by 

^i^- lW ' A^ = -4 + A;(A; + 2), (31) 

where ip = 2arctan(e :E ) and C\ are the Gegenbauer polynomials [4] (for example Cq = 1 



C\ = cosip). Thus the solution hi has one unstable mode £l = — V— with the eigenvalue 



o — ± - 

T^rif n ^nci Diffc 

cosh a; 



Actually all solutions h n (for n > 1) have an unstable mode with the eigenvalue — 1. 
To see this consider a perturbation induced by the conformal Killing vector field on S 3 
K = sinip-J^r. In terms of x-coordinate K is the generator of translations K = J^ (in fact 
this is the defining property of the coordinate x), hence 

iconf = £xh = h'(x) . (32) 

It is easy to check by differentiating Eq.(8) that £ con / satisfies Eq.(30) with X con f — — 1. 
Since by construction ^ con / has n — 1 nodes, we conclude that h n has at least n negative 



9 



eigenvalues [5]. Numerics shows that X con f = — 1 is the greatest negative eigenvalue for 
each h n and therefore h n has exactly n unstable (S'0(3)-symmetric) modes (see Table 5) []. 
This is a very interesting property which strongly suggests that the existence of solutions 

Table 5: Spectrum of perturbations 



n 


K 


K 


At 


A 4 


A 5 


1 


-i 


4 


11 


20 


31 


2 


-2.89259 


-1 


4.59106 


12.5133 


22.7280 


3 


-10.6547 


-10.6108 


-1 


4.41078 


12.6498 


4 


-94.7578 


-94.7578 


-3.06867 


-1 


4.62470 


5 


-1054.57 


-1054.57 


-10.6993 


-10.6506 


-1 



h n is due to some topological Morse-theory mechanism. This problem is currently under 
investigation. 

It is of interest to study how the energy behaves under the action of the conformal 
transformation. Given any solution h(x) a conformally deformed configuration h a is given 

by 

h a (x) = h(x + a). (33) 

For the identity solution one has h™ = — 7r/2 + 2arctan(e x+a ). This one-parameter family 
of conformally deformed configuration defines a path in the function space which passess 
through hi for a = and goes (non-uniformly) to the vacua h = ±7r/2 for a — > ±oo. It 
is easy to compute the energy along this path 

12tt 2 



1 + cosh a ' 
which, of course, for small a reproduces the linearized stability result 



£? ~ E x 



a 



dx 



cosh x 



(ti 



1\2 



6tH 



a 



(34) 



(35) 



1 It might seem from Table 5 that some eigenvalues are degenerate, in fact they form exponentially close 
pairs, which is a characteristic property of a double-well potential with high barrier. I thank H. Grosse 
and T. Hoffman- Ostenhof for pointing this out to me. 
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Note that the path of conformally deformed configurations is the path of steepest descent 
of energy from the saddle at h±. 

I have not studied non-symmetric (i.e., depending on angles ($, </?)) perturbations. 
However for the identity map hi it is known that it has exactly four unstable degenerate 
eigenmodes with A = —1, which are induced by four linearly independent conformal 
Killing vector fields on S 3 [6,7]. 

5 Proof of existence 

The detailed proof of existence of harmonic maps described above will be published 
separately. In this section I shall only sketch the main idea of the proof without going 
into the technical details. 

As I already mentioned, the asymptotic behaviour of solutions is determined by the 
autonomous system 

h" - ti + sin(2/i) = . (36) 

This equation has no regular solutions, that is, an orbit which runs into one critical saddle- 
point, say (h = — 7r/2, h' — 0), when propagated backwards, cannot reach another critical 
saddle-point {h = n/2, h! = 0) because it gets trapped by a focus at {h — 0, h! — 0) and 
spirals infinitely many times around it. It is useful to regard the coefficient tanh(x) in 
front of h! in Eq. (8) as a regularizing term which doesn't change the behaviour of orbits 
of Eq.(36) for large x but regulates the behaviour at x = by changing the sign of the 
friction term at x — 0. Actually the observation that for large x Eq.(8) is a perturbation 
of Eq.(36) might be a basis of the existence proof (cf.[8]). However this would require 
making reference to some strong results on structural stability of phase flow which I 
prefer not to invoke for a good reason that the problem may be tackled in much more 
straightforward manner. 

I start with elementary a priori properties of regular solutions of Eq.(8). 

Lemma 1 There exist two one-parameter families of local solutions of Eq. (8) near x = 
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analytic in x, b, d such that 

h( x ) = bx + 0(x 3 ) , (37) 

h(x) = d + 0(x 2 ), (38) 

where (37) and (38) correspond to odd and even solutions respectively. 

Proof. Since x = is a regular point of Eq.(8), this is a standard textbook theorem. 
Hereafter the solutions starting with initial data (37) or (38) will be called 6-orbits or 
o?-orbits, respectively. 

Lemma 2 There exits a one-parameter family of local solutions of Eq.(8) near x = oo 
analytic in c and 1/x such that 

±h(x) = -- + ce- x + 0{e- 3x ). (39) 

Proof. This fact is nontrivial since x = oo is a singular point of Eq.(8). It may be proven 
using the desingularization technique described in Proposition 1 in [8] . A b or rf-orbit will 
be called a connecting orbit if it is regular for all x > and satisfies (39) for x — »■ oo. 

Lemma 3 b-orbits and d-orbits exist for all x > 0. 

Proof. Integrating Eq.(23) from zero to some x > yields 

I r x dr, 

h'(x) = ti{0) - / — =— sin(2/i) . (40) 



cosh x Jo cosh x 

Thus h'(x) is finite for any x > and thereby h(x) is also finite. 

Lemma 4 If h(x) leaves the strip (— tc/2, it/2), then it monotonically tends to ±oo. 

Proof. Define a function 



Using Eq.(8) one has 



W(h,ti) = -ti 2 + sin 2 h. (41) 



^ = ttmhxti 2 , (42) 

dx 
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hence W increases with x > 0. Suppose that h(x ) = n/2 and h'(x ) = a > for some 
x > 0. Then for x > x 

W{x) - W(x ) = -h' 2 {x) + sin 2 h{x) - -a 2 - 1 > , 

which implies that h' 2 (x) > a 2 . 

Corollary 1 There are no regular solutions of Eq. (8) with homotopy degree bigger than 
one. 

Lemma 5 // h(x) G (— 7r/2, 7r/2) /or all x > 0, and has a finite number of zeros, then 

h! — > and /i — > ±7r/2 as x — > oo. 

Proof. First note that within the strip (— 7r/2, 7r/2) /i(x) has no positive minima nor 
negative maxima, as follows immediately from Eq.(8). Thus for sufficiently large x the 
solution /i(x) is monotonic and therefore has a limit as x — ► oo. This implies that h'(oo) = 
0, hence by Eq.(8) /i(oo) = ±7r/2 or 0. To complete the proof one has to show that the 
case h(oo) = is impossible. To prove this suppose that h(oo) = and xq is the last 
extremum of h. Without loss of generality one may assume that h has a maximum at 
Xq, hence h(x) > for all x > Xo- Now, multiply Eq.(23) by cosh and integrate by parts 
from xq to infinity. This gives 

/ — -— sin h (h' 2 + 2 cos 2 h) = , (43) 

Jxo cosh a; 

which is a contradiction unless h = 0. 

Remark. In fact Lemma 5 is true without assuming that h has a finite number of 
zeros, however then the proof is more complicated (and this stronger result is not needed 
for the proof of existence) . 

Now, having established the elementary properties of solutions, one can proceed with 
the proof. The main idea is taken from the proof of existence of globally reglar and black 
hole solutions of Einstein- Yang-Mills equations given recently by Smoller, Wasserman 
and Yau [9,10]. 
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For concreteness I shall outline the proof in the case of odd solutions (the proof of 
existence of even solutions is completely analogous). After [9] I introduce the following 
notation. Define the region T by 

r = {(h, ti, x) : \h\ < |, x > 0, (h, ti) ^ (0, 0)} . (44) 

Next, let x e (b) be the smallest x > at which the 6-orbit exits I\ For any 6-orbit define 

9(x, b) by 

0(0, b) = -- and 9(x, b) = arctan ( , fj ) forrr > . (45) 

2 y h(x, b) J 

The rotation number, Q(b), of the 5-orbit is given by 

n(b) = -±(e(x e (b),*>) + l). (46) 

Note that for connecting orbits Q = n/2 with n odd. 

Theorem 1 For each odd n there exists a connecting b n -orbit with Q(b n ) = n/2. 

The proof is based on three technical propositions. 

Proposition llfb> \pl, the b-orbit must exit T through h = ir/2 with Q(b) < 1/2. 

Proof. For b > v2, W(0) > 1 and therefore W(x) > 1 for x > which implies that h'(x) 
is positive for all x > 0. 

The following two technical propositions are crucial for the argument. They are ana- 
logues of Proposition 3.5 in [10] and Proposition 3.4 in [9], and their proofs follow closely, 
with gross simplifications, the proofs of these propositions. 

Proposition 2 Given any N > 0, there is an e > 0, such that if < b < e, then 
Sl(b) > N. 

This is a compactness property which guarantees that 6-orbits with bounded rotation 
have b greater than some positive constant. 

Proposition 3 Suppose that a b n -orbit is a connecting orbit with il(b n ) = n/2. Then for 
sufficiently small e, for b n — e < b < b n , a b-orbit exits F transversally through \h\ = -n/2 
with n/2 < n(b) < l + n/2. 
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This says that orbits which are close to connecting orbits (with b < b n ) exit T without 
making another rotation. 

Now, we are ready to complete the proof of Theorem 1. The main idea is to construct 
connecting orbits step by step. 

Step 1: Define 

7T 1 

b\ = inf {b : 6-orbit exits T via h = — with Q(b) < -} . (47) 

Proposition 2 guarantees that b\ > 0. Now, the 6i-orbit cannot exit T via h = 7r/2 because 
the same would be true for nearby orbits with b < bi, violating the definition of b\. Thus 
the 6i-orbit stays in V for all x > 0, hence by Lemma 5 it is a connecting orbit and 
fi(6i) = 1/2. 
Step 2: Let 

bi = min {b\ : &i-orbit is a connecting orbit with 0(6i) = -} . (48) 

Remark. This step is necessary because I was not able to show the uniqueness of the 
&i-orbit. 

Take b slightly smaller than b\. By Proposition 3 the 6-orbit exits V through h = —tt/2 
with 1/2 < n(b) < 3/2. Define 

71 3 

^3 — i n f {b '■ fe-orbit exits T via h — with Q(b) < -} . (49) 

By the completely analogous argument as above one concludes that 63 defines a connecting 
orbit with 0(63) = 3/2. The subsequent connecting orbits are obtained by the repetition 
of the above construction. 
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Figure captions 

Fig.l Odd solutions for n — 1,3, 5, 7, 9. 
Fig.l Even solutions for n — 2, 4, 6, 8, 10. 



